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■ We establish the exact renormalization group equation for the potential of a 

• one quantum particle system at finite and zero temperature. As an example 
we use it to compute the ground state energy of the anharmonic oscilla- 

^ I tor. We comment on an improvement of the Feynman Kleinert's variational 

^ ' method by the renormalization group. 
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Abstract 



1 Introduction 

Since most of the path integrals cannot be computed exactly, different meth- 
ods of approximation have been developed. The perturbation expansion is 
the most familiar but it diverges for all couplings strength. An another use- 
ful procedure found by Kleinert is the convergent variational perturbation 
expansion. In field theory the 'exact' Renormalization Group (RG) equations 
like the Wegner and Houghton one [0 have led to numerous non-perturbative 
results. In the path integral approach a one particle system is similar to a 
one dimensional spin lattice H, then the RG method could be also useful to 
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extract some non-perturbative results for this system. In this paper we de- 
rive the RG equation for the potential of a quantum particle which at finite 
temperature is not a closed expression contrary the field theory case. We 
show also that we can improve the Feynman-Kleinert's variational method, 
but cannot compete with the efficiency of Kleinert's systematic variational 
perturbation theory. 



2 Renormalization Group Equations in Quan- 
tum Mechanics 

Consider the euclidean action of a quantum particle at a finite temperature 



with M the mass, V the potential and T — h(5. 

The effective potential is defined as a constrained path integration over pe- 
riodic paths with period T — hp 

exp i-pVoixo)) = J Vx5{x - xo) exp (^-ls{x)^ . (2) 



where x — J dtx{t). 
The partition function is 



M 



We consider the Feynman path integral with a discretized time tn = = ne 
with an arbitrary large number, and n = 0, . . . , + 1. The Fourier de- 
composition of a periodic path x{tn) contains only a finite number of Fourier 
modes 

1 ' 

where Yl' is from 1 to y if A?" is even and from 1 to if N is odd. The x^ 

2 2coc ^"^^ 

are the Fourier modes and o;^ = 2^^^. 
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The discretized partition function is then 



dXo f-^fdXjndXm f 1„ 

/ \1{ — — ^ — exp -Y^N 

2nhe J 221^ \ h 



where the discretized action is given by 



f N+l 



n=l 



Using the fact that U^^ e^cj^ = VN + 1 and T = {N + l)e, we define the 
effective potential as (see |1]]): 

exp i-pVoixo)) = f nf ^^^^ exp (-Is^ . (4) 

Instead of evaluating (H) in one step we integrate mode after mode. Let us 
denote Vm the 'running potential' obtained after ^ — m integrations, and Vn 
the initial potential V . To find the potential Kn-i with respect to Kn we 
consider paths with only one Fourier mode: 

x{in) = a;o + exp{iujmtn)xm + H.C. 

y I\ -\- L 

We thus obtain the relation: 

exp {-PVm-l{Xo)) 

__^exp -- M^„|x„| +E^"^(^o+ , +H.C. 



(5) 



Expanding the potential Vm around the point Xq and summing over n, yields 
the following RG equation for the potential 



kipi+...+knp„=n 



m y 



(6) 
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where all the Pi >2 are different and the propagator Pm = • 

Note that in this derivation we have ignored the evolution of higher 
derivative interactions. The equation has not the closed form of the 
'exact' WH differential equation in field theory which takes into ac- 
count only the one loop contributions resumed in a logarithm term. In 
field theory the higher order loop terms are negligible even at finite tem- 
perature m but their influence is taken into account through the running 
of the parameters of the theory. In the one loop contributions are re- 
sumed in the first logarithm term but the second one contains the higher 
loop terms. To compare with the usual perturbation expansion we can 
rewrite the RG equation in the familiar cumulant expansion. Let us de- 
fine Am = (^En=o (Kn(a;(tn)) " Vm{xo) - lV^\xo)x'^{tn)) , and dcuotc the 
expectation value with respect to the gaussian weight Mtu"^ + Vj;^\xo) by 
< >m- It is straightforward to derive: 



i-ki 



E/ 1 \fcl + ... + fcn V \ ^/ J \ m \ ^/J 

kjpi+...+{l-k-,...-k„-i)p„=n l^l-yPl-) l^n-{Pn-) 

Then defining the cumulants by 

< > =< (A - < A > r > 

the Renormalization Group equation becomes: 

Vm-l{x,) = V;,(xo) + ^logfl+^^'^^''° 




(3 "V ^IM 

In this formulation the equation has the disadvantage to mix the powers of 
1 

Clearly the RG equation seems not to be useful to extract nonperturbative 
results. Actually it contains an infinite number of terms, negligible only if 
are small (for finite temperature the coupling constants must be small). 
Nevertheless it is interesting to discuss the various limiting cases. 



4 



2.1 Infinite temperature 

For a large temperature the quantum fluctuations are small and the system is 
close to the classical one. The flow of the potential stays near the classical one 
(the potential energy of the action) and will be obtained after a relatively 
small number of iterations of the RG equation. In particular in the limit 
^ _^ m.^ 1^0) g^j^j^ p g^j^g q£ order and we obtain 



a running potential constant along the flow. As a consequence the quantum 
partition function reduces to the classical one. 

2.2 Zero temperature 

In this opposite case due to quantum fluctuations the effective potential is 
expected to be different from the classical one and will be obtained after a 
huge number of iterations of the RG equation (typically 10® for (3 = 10^). The 
higher loop contributions are negligible, then we can extract informations for 
large coupling constants. In particular, in the limit (3 ^ oo the propagator 
Pm is of order —(^y — . Then keeping only the one loop contributions in (|^) 

we obtain the closed form: 



Recall that for large A^, cu^ is equivalent to ( ) and then runs from to 



Vm-l{Xo) = Vm{Xo) 








the limit N ^ oo gives a continuous equation 
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where we have introduced the notations k"^ = tu^, Ak = ^^f^, and neglected 
the higher loop contributions which are of order (Afc)^. This equation is the 
one dimensional Wegner Houghton equation for the effective potential [0]. 
It is useful to extract informations even for very large coupling constants as 
shown below. 



3 Examples 

We solve exactly the Renormalization Group equation for some easy examples 
and recover known results. 



3.1 The free particle 

This case is trivial because Viv(a;o) = 0. It's easy to show recursively that at 
each step Kn_i(xo) = Vm{xo) = 0. Then, the effective potential Vo(a;o) = 0. 

3.2 The harmonic oscillator 

As a second example, consider the case: Vn_{xo) = y^'^^o- O^ice again, we 
obtain recursively: Vm-i{xo) = ^ii^x^ + Cm with Cm a constant. Inserting 
this result in (H), we obtain the solution 

M . 2 ,1, f^K (ul + Q?' 
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V,{x^) = —n'xi + -\og\Ili 

Z, p y y UJ, 

which gives the partition function (see ]l|): 

1 ( \ 

" 2 I m \ 



3.3 Perturbation expansion 

Replacing Vm by Vn in ([7p one recovers the usual perturbative expansion 



Voixo) = V.{xo)+Y.l^'^og\l + ^—-\+j-Y.<A.>. 



2 I I Vn (xn) \ 1 2 



" iulM hp 



m=l 
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JV 

4 The anharmonic oscillator 

In this part we compute the ground state energy of the anharmonic oscillator 
at zero temperature for various values of the coupling constant by using the 
RG equation, and compare with Kleinert's [|l|] variational results. The ground 
state energy is the value of the effective potential Vb(xo) at its minimum 
xq = 0. Unfortunately it is a hard task to use the RG equation (||) because 
at each step we must find the second derivative of the potential by fitting it 
numerically. Due to the very slow convergence we must repeat the procedure 
a huge number of time for each point of the potential. We leave this method 
for an another paper and try instead to compute the flow equations for the 
coupling constants. 

We define the n-th coupling constant at the scale m by gl^^ = -f^ Vm(xo). 

"^0 1x0=0 

We give a general formulation for the flow of the couplings where we neglect 
terms of higher order in ^ : 

9^:u=9i^^ + bog(i+^]. (10) 



For A; 7^ we have 



\ P \ai+...ap=k,a,>0 ^ "l' "P" 

(11) 

{ai+2) (ap + 2) 

where Ap is the combinatorial factor of ^^^^ — . . . ^^^i — in the series expan- 
sion of (V^'^\xo)Y in powers of Xq and Pm(0) = — — - — m- 

m 

For the initial potential we choose 

Vn (x) = -Un^x" + (12) 
2 ^ ' 2 4! 
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9m=0 


ym=o 


i/m=0 


0^'^ 
i/m=0 


ym=o 


Erg 


1.5140 


2.4 











0.56134 


1.4522 


1.5343 











0.55807 


1.4643 


1.7706 


3.9484 








0.55855 


1.4617 


1.7110 


2.8109 


-16.662 





0.55847 


1.4662 


1.7236 


3.0112 


-16.794 


-237.46 


0.55848 



Table 1: comparison of the ground state energy Egr of the anharmonic 
oscillator for M = Q = 1, X = 2.4 and different truncations of the potential. 
The exact value is E exact = 0.55915 



The numerical solution of the flow equations (|T0|) (0) shows a good conver- 
gence for = 10^ and (3 = 10^. There are no significant improvement for 
larger values of iV or /3. At first, we have tried different truncations of the 
polynomial potential and found the best values by keeping only the couplings 
g^'^^ and g^^^ as is shown in table 1 where we can also read the different values 
of the couplings at the end of the flow when m = 0. Adding more coupling 
does not improve the result, on the contrary. In fact in order to get better 
results we can not truncate the effective potential because in quantum me- 
chanics each coupling constant is relevant. Even when the couplings g^^"^ and 
g^^^^ take large negative values they do not influence very much the flows of 
the other couplings. Then, it seems that around its minimum the potential 
is well fitted by a polynomial interaction of order six. 

A detailed comparison of the truncation at with Kleinert's variational 
method is given in table 2. All the couplings are relevant because they grow 
as we iterate the RG equation. Surprisingly, the truncation of the potential 
at the six order gives always better results than those of the variational 
method. The advantage of the latter lies on the fact that one can improve it 
systematically with a perturbation expansion which is convergent (see ||l[). 
In our case the only way to improve it is to consider the potential as a whole. 
This requires a numerical study which is difficult to achieve due to the very 
slow convergence of the iterations in the zero temperature case. 
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A 


Eexact 


E 


Erg 


(2) 
9m=0 


(4) 
9m=0 


J6) 
9m=0 


2.4 


0.55915 


0.5603 


0.5585 


1.4643 


1.7706 


3.9484 


4.8 


0.60240 


0.6049 


0.6014 


1.8189 


3.2687 


12.078 


7.2 


0.63799 


0.6416 


0.6366 


2.2126 


4.6951 


22.058 


9.6 


0.66877 


0.6734 


0.6672 


2.4041 


6.0842 


33.296 


12 


0.70618 


0.7017 


0.6943 


2.6614 


7.450 


45.527 


14.4 


0.72104 


0.7273 


0.7190 


2.9031 


8.800 


58.589 


16.8 


0.74390 


0.7509 


0.7417 


3.1324 


10.137 


72.384 


19.2 


0.76514 


0.7721 


0.7628 


3.3514 


11.466 


86.832 


21.6 


0.78503 


0.7932 


0.7825 


3.5619 


12.787 


101.87 


24 


0.80377 


0.8125 


0.8011 


3.7648 


14.102 


117.46 


240 


1.50497 


1.5313 


1.4982 


14.735 


128.19 


2525 


1200 


2.49971 


2.5476 


2.4877 


41.683 


627.76 


21462 


2400 


3.13138 


3.1924 


3.1162 


65.742 


1250.3 


54006 


12000 


5.31989 


5.4258 


5.2937 


190.83 


6222.1 


460992 


24000 


6.69422 


6.8279 


6.6611 


302.50 


12433 


1161228 



Table 2: comparison of the exact ground state energy Eexact with the varia- 
tional energy E^ar and the energy Erg obtained from the flow of our trun- 
cated potential. 
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A 


Erg 


m 

9m=0 


ym=o 


ym=o 


4.8 


0.073 


0.366 


0.821 


4.288 


7.2 


0.189 


0.681 


1.955 


13.02 


9.6 


0.266 


0.960 


3.111 


23.34 



Table 3: Ground state energy and flow of the coupling constants for the 
double well potential. 

4.1 The double well potential 

In this part we consider the case of the double potential with a frequency 
Q = —1. The results for some values of the coupling A (see table 3) show that 

(2) 

grnLo > 0. It means that the quantum fluctuations smear out completely the 
double well as it must be, because the effective potential is a convex quantity 
at T = 0. For A = 2.4 our truncation is not accurate enough to get a 
positive term in the logarithm of equation (H). For all larger couplings our 
approximation is enough to show the spreading of the double potential by 
the quantum fluctuations. 

5 The Feynman Kleinert method 

In this section, we give some comments concerning the FK method improved 
by the renormalization group. In the FK method, one tries to flnd a quadratic 
potential at each point Xq fitting at best the effective potential. One im- 
proves this procedure by looking for a quadratic potential Mflrn^^ fitting the 
potential Vm{xo) at each step of the renormalization group fiow. Then by 
improving in such manner the FK method one will take into account some 
contributions of the Kleinert's variational perturbation expansion [|[]. To 
achieve this, we follow Kleinert's method and insert in the right hand side of 
(ID the trial frequency flf^{xo) (see |jl[]). 
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e 



Then, using the Jensen-Peierls inequahty 

c?/i(x) exp (— 0(x)) > exp ( — / dfi{x)0{x] 



for any positive measure /i normahzed to one, and following the same steps as 
in Kleinert's book we obtain the following RG inequality for the potential 

Vm-l{Xo) < 

where 

2r ^ 2 1 



Let us call \42^ the integral on the right hand side of the inequality. Minimiz- 
ing the right hand side in the variable Vt^ we obtain the best approximation 
for Vm-i- The frequency 

gives the following equation, which can be solved recursively, starting from 

= 

nlixo) = /^exp (-xV2) -^Vi'\xo + ^al{x,)x). (14) 

J \/27r ^ ^ 0. ^ 



"271 ^ ^ 2^a. 

Keeping this value for fi^, enables us to rewrite the approximation 

V^..M = ^log ^1 + — j - -^^^TJ^^^^ + (15) 

which is the variational RG equation. For finite f3, it is possible to solve 
(|15D, firstly by computing recursively f2^(xo), secondly by calculating V^2^ 
and inserting this result in the Renormalization group equation. We plan to 
use this equation in an another paper. In the following section we check its 
validity in the limiting cases. 
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5.1 Infinite temperature 

Recall that for large N, u;^ is equivalent to (^)^ Then for (3^0, 

and „2^i" 2 are of order and aL(xo) is of order /3. Inserting these results 
in the equation for V^2^ gives V^2^ = V^„(a^o)- For /? — > we recover the usual 
invariance for the flow of the effective potential: 

Vm-l{xo) = Vm{xo) . 



5.2 Zero temperature 

In the limit /3 ^ oo, a'^{xo) is now of order ^. To obtain we expand the 

M 



potential Vm in a series of a^. To lowest order we obtain: il'^{xo) = ^"^ 



In the same way, we get the lowest order: 1^2^ (xq) = Kn(2^o)- 
Then keeping only the relevant terms, the limit for /3 — > gives again the 
Wegner Houghton equation: 

V„-,W = V„(.o)4log(l + ^) 

To understand why we recover the WH equation, one compares with the 
Kleinert's variational perturbation expansion. In the RG equation for the 
potential at T = 0, the contributions of all the higher loops terms are included 
in the flow of the couplings. By choosing a variational frequency which flows 
in the RG equation we automatically take into account all the higher loops of 
the variational perturbation expansion and the result must be independent 
of the variational frequency. In other words, by integrating mode after mode 
we automatically resume the perturbation expansion and we cannot improve 
it. A possibility to resume partially the higher loops would be to integrate 
a lot of modes at each step, and to deduce the variational frequency for the 
resulting potential. 



6 Conclusion 

The RG equation for the potential of a quantum particle at finite tem- 
perature was derived. This equation does not allow us to compute non- 
perturbative quantities. Then we suggest to use the Feynman Kleinert's 
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variational method improved by the RG. At zero temperature we recovered 
the Wegner Houghton equation which was used to compute the ground state 
energy of the anharmonic oscillator. It would be also interesting to compute 
the RG equation for the effective action. Our preliminary work show the gen- 
eration of non local interactions. We also plan to apply the RG procedure 
for systems with many quantum particles in the discretized path integral 
representation. 
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